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Ultimately,  et  relatively  large  dietanoes  from  the  eouroe,  there  ia  some 
similarity  between  the  blast  waves  from  all  explosions.  Various  authors  have 
demonstrated  this  mathematioally,  but  have  not  drawn  attention  to  the  feet  that 
their  results  are  not  entirely  compatible  with  the  olaaaieal  hypothesis  that, 
ultimately,  blast  waves  become  less  and  less  aware  of  the  nature  of  the 
explosion  and  tend  to  depend  only  on  the  magnitude  of  the  energy  release. 

The  same  hypothesis  has  been  used  to  derive  atmospherie  scaling  laws,  whloh 
therefore  need  alternative  interpretation  if  the  hypothetic  ia  found  to  be 
falee. 


The  ultimate  flow  behind  a  decaying  spherical  shock  ia  studied  here  from 
this  point  of  view.  It  is  conclule-l  that  the  classical  hypothesis  is  false, 
and  la  not  aasentlal  to  the  matheoatioal  similarity.  It  has  been  shown 
elsewhere  -tSoJ  that  It  it  also  not  essential  to  the  atmospheric  scaling  rules. 


The  theory  developed  in  making  these  conclusions  suggests  that  the  flow  < 

behind  a  spherical  shook  ia  divided  into  two  parte  by  an  inner  sphere,  such  \ 

that  the  flow  within  this  sphere  has  no  influence  on  the  deosy  of  the  shook. 

The  excess  energy  in  the  outer  part  increases  indefinitely  In  proportion  to  j 

the  shook  radius,  and  that  within  the  inner  sphere  must  therefore  correspondingly 
decrease  indefinitely  since  the  total  excess  energy  oust  remain  finite  and  { 

equal  to  the  energy  released  by  the  explosion. 

A  hypothesis  is  suggested,  in  an  attempt  to  obtain  a  higher  order 
asymptotic  solution,  for  the  blast  wave,  than  the  well-known  form  given  by  \ 

previous  authors.  This  solution  is  coopered  fully  with  experimental 
observations.  i 
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1.  Introduction 

It  has  long  been  known  tl<at,  far  away  from  tha  origin  of  explosion,  there 
la  ultimately  none  similarity  between  the  bloat  waves  from  all  explosions, 
whatever  the  source  and  nature  of  the  energy  release  which  originally  Initiated 
them.  Be  the  [ 1 ]  first  showed  that ,  if  R  is  the  radius  of  a  decaying 

apherloal  shook,  then,  as  R  "  «.,  tha  peak-overpressure  behind  tha  shook 
ultimately  behaves  aa  [.1  log^  R]  “ tha  duration  of  poaitiva  overpressure 
behaves  aa  log1'  R,  and  the  impulse  per  unit  are*  in  the  poaitiva  overpressure 
phase  behaves  aa  R~  Kirlcwood  and  Brinkley  [2]  developed  a  theory  oovering 
tha  whole  course  of  an  explosion,  baaed  on  a  physical  assumption  about  energy, 
which  led  to  the  same  asymptotic  form  for  large  value*  of  R,  later, 

Whitham  [3]  re -derived  Beths' s  results  by  an  alternative  approach,  aimed,  Ilka 
Be the' a  method,  at  correcting  tha  imperfections  of  tha  'aooustio'  solution  of 
the  problem  of  a  decaying  shock. 

Both  Bet he  and  Kirkwood  and  Brinkley  define  the  energy  of  the  shook  wavs, 
aa  it  crosses  t)te  sphe.-e  of  radius  R,  aa  the  nett  work  done  on  the 
undisturbed  atmosphere  exterior  to  the  sphere  R,  But  only  by  making 
assumptions  about  tha  motion  In  the  negative  phase  of  the  blast  wave,  and  thua 
assessing  effectively  tha  final  nett  amount  of  rcrk  done  on  the  undisturbed 
atmosphere  external  to  the  sphere  of  radius  R,  do  they.reaoh  the  conclusion 
that  this  nett  energy  tends  ultimately  to  Eero  as  log's  R  when  R  Is  large, 
thua  indicating  a  very  slow  final  dissipation  of  energy,  a  result  first  derived 
. ly  Penney  [11].  None  of  theoe  authors  has,  however,  drawn  attention  to  the 
difficulties  associated  with  considerations  of  tha  lraediate  energy 
distrilxition  behind  a  decaying  apherloal  shook.  If,  instead  of  considering 
the  total  nett  work  done  on  tha  undisturbed  atao sphere, external  to  tha  sphere 
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Sa  the  development  of  thiu  energy  with  tine  la  considered,  It  la  found,  at 
once,  without  any  assumptions  concerning  the  negative  phase,  that  the  work 
done,  up  to  the  end  of  the  positive  overpressure  phase,  on  the  atmosphere 
originally  outside  the  sphere  R,  behaves  as  R  when  R  is  large,  and  this 
Increases  Indefinitely  as  «  increases.  This  is  a  curious  result,  and, 
despite  the  faot  that  suoh  a  curious  energy  situation  nuat  be  restored  ty  a 
negative  phase,  it  is  of  the  greatest  Importance  to  an  understanding  of 
explosion  phenomena  (particularly  as  applied  to  target  response)  to  examine 
the  consequences  of  such  an  anomalous  energy  distribution,  Whltham  [3]  does 
not  deal  with  energy  oonsideratiors,  but  in  a  later  paper  [12]  he  mentions  a 
similar  anomaly,  namely  that  the  nasa  flowing  from  inside  the  sphere  R  to 
outside  the  sphere  R  during  the  positive  overpressure  phase  ultimately  behaves 
as  R,  and  thus  necessitates  the  existence  of  a  negative  phase  in  which  it  can 
return.  The  anomaly  does  not  exist  in  plans  one -dimensional  flow;  In  this 
oase  the  excess  energy  in  the  positive  overpressure  phase  tends  to  a  finite 
quantity  aa  R  increases.  In  axially  symetrlo  flow,  the  excess  energy  In  the 
positive  overpressure  phase  behaves  ultimately  aa  R“  when  R  is  large. 

In  the  absence  of  aiy  real  attempt  to  faoe  up  to  the  difficulties  of  a 
positive  overpressure  phase  containing  an  indefinitely  increasing  amount  of 
energy,  continuing  appeal  has  tacitly  been  made  to  the  o lassie al  hypothesis  of 
explosion  theory,  namely  that,  for  large  R,  the  motion  behind  a  blast  shock 
becomes  dependent  dominantly  on  a  finite  quantity  of  energy  associated  with 
the  energy  released  ty  the  explosion,  and  less  and  less  on  the  nature  aril  mode 
of  the  energy  release.  The  asymptotic  results  described  above  would  only  be 
compatible  with  suoh  a  hypothesis  if  the  energy  in  the  positive  overpressure 
phase  tended  to  a  finite  non-sero  amount  for  large  R,  since,  aa  will  be  shown 
In  this  paper,  the  decay  of  a  spherical  shook  is  determined  entirely  ty  the 
motion  in  a  limited  region  of  flow  just  behind  the  shod:,  not  such  different 
from  the  positive  overpressure  phase.  The  hypothesis  is  therefore  now 
believed  to  be  false. 

This  same  classical  hypothesis  has  been  used  ty  Sachs  [13]  to  derive 
atmospheric  scaling  laws,  though  the  seme  sealing  laws  were  derived 
previously  by  Taylor  [14],  olosa  in  to  a  very  intense  explosion  [of.  10]  on  the 
different  hypothesis  that  the  total  excess  energy  within  the  spherical  shook 
is  conserved.  The  present  author  [10]  has  shown  that  Sachs'  rule  far 
atmospheric  coaling  can  be  derived,  and  generalised,  without  appeal  to  the 
classical  energy  hypothesis ,  and  that  Its  practical  usefulness  ia  not  therefore 
nullified  by  the  falsity  of  the  hypothesis.  Instead,  Sadis'  rule  Is  seen  to 
be  not  an  asymptotic  law,  as  first  derived,  but  a  rule  for  the  'middle- 
distanoe'  which  ia  of  far  greater  importance  in  practical  applications. 

Taylor's  derivation  of  the  atmospherlo  sealing  rule,  olose  in  to  a  very  intense 
explosion,  also  new';.-  further  consideration,  since  the  high  temperatures 
associated  with  the  high  energy  release  per  unit  volixas  of  a  very  Intense 
explosion  imply  a  high  rate  of  loss  of  energy  from  the  system  within  the 
spherioal  shock  ty  thermal  radiation  [of.  15].  The  derived  atmospherlo 
scaling  rules  may  therefor©  only  be  expected  to  apply  approximately,  under  the 
oonditions  of  Taylor's  similarity  solution,  over  ranges  of  R  in  which  the 
energy  lost  to  the  system  ty  thermal  radiation  is  eitlier  small  eoepared  with 
the  total  energy  of  the  system,  or  scales  consistently  with  the  derived 
sealing  rules. 

In  the  present  paper,  the  equations  of  unsteady  non-homentropio 
oomprasaible  flow  with  spherioal  gymetry  are  reduced  to  oharaoteristio  form 
(sect ion  2) ,  in  terms  of  variables  whioh  are  effectively  the  Rieraann  functions 
for  the  corresponding  plane  flow.  The  boundary  oonditions  at  s  general 
spherical  shook  front  ere  derived  in  the  same  variables  (section  3)  and  an 
iterative  solution  is  developed.  This  method  of  solution  was  developed 
independently,  lut  was  found  later  to  be  similar  to  that  used  ty  Bethe  [  l] , 
although  the  present  solution  is  carried  further  than  Bethe 's  solution. 

Energy  and  impulse  considerations  are  introduced  (aection  4) ,  to  derive  the 
limitations  on  the  form  of  the  shock  decay,  anl  the  results  derived  are  in 
acoord  with  previous  results  described  above,  but  are  taken  to  higher  order 
approximations.  Tl-e  existence  of  a  limiting  characteristic  or  wave-front, 
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behind  the  leading  shook,  Is  demonstrated  ( Motion  5)»  s»  was  also  done 
previously  by  iThithsm,  namely  the  first  positive  oharaoterietio  which  fails  to 
overtake  the  shook  in  finite  dlatanoe.  This  marks  the  re«r  of  what  is  defined 
here  as  the  shook  decay  phase  of  the  motion,  for  it  divides  the  flow  behind  the 
shook  into  two  parts,  a  spherical  annulus  (between  the  shook  and  the  limiting 
characteristic)  whioh  alone  controls  the  decay  of  the  shook,  and,  inside  it,  a 
spherical  region  whioh  cannot  influence  the  shook  deoay  process.  The  previous 
results  (of  section  k)  then  show  that,  although  the  total  axoesa  energy  within 
the  leading  shook  is  always  finite  and  equal  to  the  energy  released  by  the 
explosion  (less  possible  radiation  loaaaa),  ita  components,  in  the  two  regions 
defined  above,  respectively  Increase  and  decrease  indefinitely  in  proportion  to 
R ,  when  R  la  large ,  and  thus  oontradiot  the  olaasloal  hypothesis  that  the  flow 
in  the  shock  decay  phase  immediately  behind  the  leading  shock  is  determined 
dominantly  by  a  finite  quantity  of  energy. 

The  difficulties  of  obtaining  a  higher  order  approximation,  without 
recourse  to  hypothesis,  or  to  the  flow  conditions  behind  the  shock  decay  phasa, 
are  introduced  (section  6),  and  a  hypothesis  for  a  higher  order  approximation 
la  suggested.  On  such  a  hypothesis,  the  ultimate  form  of  the  pressure-pulM 
In  a  blast  wave  ia  derived  ( section  7) ;  the  higher  order  results  from  the 
hypothesis  are  compered  with  experimental  observations  (Motion  8) ,  and  show 
promising  agreement.  Nevertheless,  even  If  these  higher  order  results  are 
Ultimately  found  to  be  sound,  they  do  not  apply,  with  practical  accuracy, 
closer  In  than  about  20  ft.  from  a  nna-pound  charge  of  conventional  explosive 
at  sea-level,  where  the  slock  peak  overpressure  is  about  one-sixth  of  an 
atmosphere.  It  still  retains  true,  therefore,  as  stated  by  Beths  [l],  that 
"there  is  considerable  danger  In  using  theee  relations  for  moderate  pressures 
where  the  pulse  has  not  yet  reaohed  ita  limiting  form". 

The  real  practical  problem  of  explosion  theory  therefore  remains,  namely 
to  provide  a  consistent  and  oastprahenaivs  working  solution  for  the  'middle- 
diatanoe'  of  moderate  peak  overpressures.  For  it  ia  in  the  ' middle-distance' 

[ of.  1 6]  that  the  oritioal  oonditlona  of  target  response  generally  ooour. 

This  paper,  then,  makes  no  contribution  to  the  really  important  practical 
problem,  except  in  so  far  as  it  Maks  to  improve  understanding  of  the  ultimate 
asymptotic  form  of  a  bleat  wave,  which  must,  of  neoesaity,  be  an  end-oondition 
for  the  'middlo-di  stance '  problem.  It  Is  dear  that  a  proper  understanding  of 
the  end-condition  ia  a  neoessary  preliminary,  In  any  ease.  It  may,  or  may 
not,  provide  the  key  also  to  a  useful  appro aoh  to  the  more  practical  problem. 

2,  Unsteady  non-homentroplo  flow  with  gpherloel  wnsetry 

The  equations  of  motion  of  the  air  behind  the  leading  bleat  shook  in  a 
spherical  explosion,  in  regions  where  no  further  shocks  ooour,  are) 
from  the  conservation  of  momentum, 

Ut  ♦  uUj.  ♦  Pj/p  ■  0}  (2.1) 

from  the  ooneervatlan  of  mats, 

pt  *  uflj.  ♦  pup  ♦  2up/r  ■  0;  (2.2) 

and  from  the  oonservstion  of  energy 

Pt  +  upp  b  a*(»t  ♦  uPr).  (2.3) 

p  derates  pressure;  p,  density;  a,  sound-spaed;  u,  partlole  velocity;  r 
derates  distance  from  the  oentre  of  spherical  syirmetry  and  t,  time;  suffixes  r 
and  t  denote  partial  derivatives;  and  air  is  assumed  to  behave  as  an  ideal 
polytropio  gas,  with  aquation  of  state 

{$■ 


and  sound -speed  a 


ar-m- 

■  (yp/e)^» 


(2A) 

(2.5) 
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S  denotes  ant ropy ,  a,  epeoifio  heat  at  eonetant  volma*;  and  y,  the  ratio  of 

epeoifio  heata  at  eonatant  pressure  and  oonatant  volute  re»p*ctlv*ly,  namaly 
Oj/oy,  la  taken  to  have  the  value  1/ 5  for  airt  suffix  o  refer*  here  to 

oonditiona  In  the  undiatuxhed  uniform  atmoaphere  at  rest.  In  addition 


*T  "  y-T^ 

(2.6) 

«a& 

f  -  <X* 

(2.7) 

where £  la  a  oonatant,  and  •  denotaa  ahaolute  taoperature. 

Fundamental  unite  are  now  ohoaen,  defined  tty  the  three  quant  it  lea,  atmospheric 
pressure  pg,  ataoapherlo  sound-speed  a0,  anl  an  arbitrary  length  L,  and  the 
following  notation  la  adopt ad, 


1,0.  In  general  £  ■  P/tP] ,  where  P  denotea  ary  phyaloal  quantity  and  [p]  la 
that  combination  of  the  fundamental  unita  which  ha*  the  eeme  dimensions  ea  P, 
In  addition,  the  oomrentlon  la  adopt'd  of  writing  3  ■  (S  -  30)/ft.  With  this 
notation  and  ualng  (2,6) ,  (2.4)  and  (2.5)  m*T  ho  written 

£  .  a”/<*  “  0  ."ft  (2.8) 

f  .  y4a/(v-r)  #-i  (2.9) 

The  equation*  of  motion,  (2.1),  (2.2)  and  (2.3),  n ay  now  he  written  In  the  new 
notation,  end  p,  £  may  he  eliminated  from  them  ty  mean*  of  (2.8)  and  (2.9). 

The  reault  la  ** 

*t  ♦  2  Jr  ♦  *  %  *  **  7^1}  %  ■  0  (2.10) 


St  +  £  Sr  *  0  (2.12) 

In  whioh  X*  -  (y  -  l)/(y  ♦  l). 

Mew  dependent  variable*  «  and  fi  are  now  ohoaen  instead  of  u  end  a,  auoh  that 

a  »  (1  -  X*)(a  -fi),  (2.13) 

a  -  1  ♦  X*  («  +  /»),  (2.14) 

or*  altematiTely 


4* 


(2.15) 

2X*  2(1  -  X*)* 

(a-1)  .  u 

2X*  2(1  -  A*)* 

(2.16) 

from  whioh  it  may  be  teen  that  a  and  p  are,  effectively,  the  Rienam  function* 
whioh  are  constant  along  eharecterlstie  lines  in  unsteady  plane  hooen tropic 
fixer,  She  equations  then  reduce,  after  suitable  manipulation,  to 

a,  ♦  [n«-  0  -sx«)j]  s  ♦  ,  0  (2.17) 

*  £  £  2(1  ♦  X*J  — 


fit  -  tl-(l-2X*)ae/»]  fir  + 


(«-/9)[1  ♦  X*(«+/ 


2(1  +  X* 


fit  ♦  (1  -  X*)(«  -Jjfig  -  0 


0  (2.18) 
(2.19) 


In  this  form  it  is  seen  that  the  system  of  equations  is  hyperbolic,  and  ha* 
three  families  of  ohareoteristio  line*  in  the  (r,  -plane.  Along  the  positive 
oharaoteristio  lines, 

d£  .  [1  ♦  «  -  (1  -  2XV)dt,  (2.20) 

the  intermediate  integral  relation  satisfied  is,  from  (2.17) , 

«. . i0i s . ,  .  o.  (,.«> 

(1  ♦  «  -  (1  -  2X*)rt  £  2(1  ♦  X*) 

Along  the  negative  charenterlatlo  lines, 

<fe  ♦  (1  -  (1  -  2X*)«  e  0]  -  0,  (2.22) 

the  intermediate  integral  relation  satisfied  is,  from  (2.18), 


j.  (a  -  *)[1  a  X»(«  ♦  /ft]  «£.  ti»X»(«ed 
’  (1  -  (l  -  2X*)«  ♦  Pi  £  2(1  ♦  X*) 


dS  a  Oj 


and  along  the  world-lines ,  or  particle  paths  in  the  (r,  t) -plane, 

dr  -  (1  -  X*)(«  -  fi)  dt,  (2.24) 

the  relation  satisfied  is,  from  (2.19), 

S  a  oonstant.  (2.25) 

J.  The  decaying  spherical  shook 

When  the  leading  blest  shock  from  a  spherioal  explosion  has  reached 
distances  large  occpared  with  the  else  of  the  explosive  oharge,  or  when  the 
peek  overpressure  lmediately  behind  it  is  email  oorapared  with  the  undisturbed 
atmospheric  pressure,  it  is  pertinent  to  ask  whether  the  actual  conditions 
under  whioh  the  energy  was  released  by  the  'detonation'  continue  to  have  ary 
influence  on  the  way  in  whioh  the  shock  decays,  or  whether,  in  feet,  the  prooesa 
of  decay  tends  to  depend  less  and  lesa  on  the  exact  nature  of  the  energy 
release.  The  answer  to  this  question  entails  the  study  of  an  lnoonpletely 
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formulated  problaa  In  whioh  the  governing  equation*  are  those  of  Motion  2 
above,  and  the  relevant  boundary  condition*  are  thoae  at  the  leading  bleat 
shock,  and  the  end-oondltlon  that  the  (hook  ultimately  deoaya  to  aero  strength 
as  Its  radio a  increases  Indefinitely.  The  problem  la  to  draw  conclusion* 
from  this  data  only,  without  recourse  to  the  Initial  conditions  of  th* 
detonation,  or  any  other  oonditiona,  at  any  tine,  at  boundaries  nearer  to  the 
eentre  of  the  explosion  than  the  leading  shock. 

Let  the  position  of  the  leading  blast  ah>ck  at  any  time  be  denoted  ly  th* 
equivalent  relation*  R  ■  R(t)  or  X  ■  T(r).  In  th*  notation  of  aaetlon  i! 
above,  the  t ine -dsr elopnsnt  of  the  leading  blast  shook  nay  b*  denoted  ty 
T  •  T(R) ,  and  it  la  oonvenisnt ,  whan  th*  leading  shook  la  weak,  to  write  this 
relation  In  th*  font 

£*&-•&)  (3.1) 

where  *(&)  la  an  unknown  function  of  g,  to  be  determined  so  far  as  la  poaeibl*. 
Than  fron  (3.1),  if  dashes  denote  differentiation  with  reapeot  ta  £, 


H 


(3.2) 


la  the  Maoh  laasber  of  the  ahoak  In  terms  of  the  speed  of  sound  In  th* 
undisturbed  etna  sphere,  and  thus,  for  a  weak  ahoak,  *'  la  email  compared  with 
unity.  Th*  shook  ooiriitlona,  giving  th*  state  of  th*  air  lanedlately  behind 
th*  blast  shook,  are  (see,  for  example,  (4),  pp.  74,  73), 


£  .  1  ♦  (1  ♦  X'JOt*  -  1) 


Y  1  ♦  X*(M*  -  1) 

Substituting  for  M,  from  equation  (3.2),  th*  fallowing  result*  are  obtained, 

£  B  1  ♦  2(1  ♦  X*)e*  ♦  3(1  ♦  V)a'#  4  0(e'*) 

a  *i  2(1  -  X*)«'  ♦  (1  -  x*)**1  4  0(e'*) 

-  (3J 

=  ■  1  ♦  2(1  -  X*)*'  ♦  (3  -  4X*)(l  -  X*)e'*  ♦  0(e'») 

awls  2X*  •*  e  X*  *'•  4  0(e'*)» 


,  fron  (2.15),  (2.16)  and  (2.4), 

1  ■  1  u 

2X*  2(1  -  X*) 


2*'  ♦  *'*  4  0(«'*) 


fi 


(a  -  1)  R 

2X*“  2(1  -  X*) 
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0(t'») 

0(*'»). 


(3.3) 


Th***  are  the  required  fora  of  the  leading  shook  boundary  conditions  la  th* 
notation  adopted  hare. 
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The  end -condition ,  that  the  shook  tends  to  aero  strength  me  the  ehook- 
redius  increase*  Indefinitely,  is  given  by, 

£  «  •  t  e'  •  0  •  (3.8) 

The  problem  then  is  to  find,  ao  far  as  1*  possible,  the  solutions  of  the 
governing  equations,  expressed  ohamoteriatioally  In  (2.20)  to  (2.25), 
satisfying  the  Incomplete  set  of  boundary  conditions  (3.5)  and  (3.6).  Such 
solutions  will  apply  to  the  sir-flow  Immediately  behind  the  leading  blast  abode, 
ao  far  as  no  further  shocks  occur  In  the  flow. 

The  method  of  solution  adopted  here  may  be  described  as  one  of  'snelytio 
iteration' ,  It  was  found  ultimately  to  have  been  proposed  and  used 
previously  to  some  extent  by  RJL.  Beths  [l)j  Widths*  [2]  has  adopted  a 
different  method  of  solution.  Both  methods  are,  in  fact,  devices  for 
Improving  on  the  gross  laperfeotlona  of  'aoouatlo'  solutions  to  the  problaa, 
which  fall  on  aeoount  of  the  divergence  of  the  ehereeterlstloa  at  Infinity  In 
the  proper  solution. 

It  is  proposed  to  seek  a  solution  only  so  far  as  sooond  order  tense  In  t* 
at  Boat.  Since  8,  at  the  leading  shook  la,  by  (3.5)*  only  of  the  order  of  #** 
and  la  constant  along  world-linos  whose  slope  la  of  order  It  follows  that  8 
la  of  the  order  s'*  at  all  points  of  Interest  and  nay  be  neglected  ao  far  as  a 
solution  to  aooond  order  terms  in  »•  la  concerned.  Slnoe,  moreover,  by  (3*5) , 
ft  is  only  of  order  #'*  at  the  leading  shook,  whilst  s  is  of  order  s',  a  first 
approximation  to  equation  (2.21)  is 

d(e  j)  ■  0. 

This  is  the  approximation  Introduced  by  Beths  [l ],  and  given 
ar  ■  oonstant  ■  2  ij 

or  «  »  2  if/z  (3.7) 

on  the  positive  characteristic  which  overtake!  the  limiting  shock  at  0  (Figure 
l).  Being  this  result,  a  first  approximation  to  the  equation  of  the  positive 
characteristic  (2.20)  Is 

giving,  by  Integration, 

r  -  t  -  2  «o  Ro  log  r  >  *o  ~  2  *o'  &  log  B0  (3.8) 

for  the  positive  oharaoteristie  which  overtakes  the  leading  ebock  at  0. 

(H.B.  All  logarithm#  ere  natural  logarithms  to  baas  e.] 

The  negative  oharaoterlatlo  (2.22)  which  intersects  the  leading  shook  at  A 
(Figure  l)  is,  to  a  first  approximation,  elsply 

1*  1  m  £a  +  £a  ■  2  5l  '  *A  (3*9) 

If  the  oharaoterlatlo  which  meets  the  leading  shook  st  0  outs  this  negative 
oharaoterlatlo  at  *,  then,  at  B,  to  the  mow  order  of  accuracy, 

2  £  *  2  Si  ♦  *i  ■  2*tS§0l'°sr«'«!0-2e»80lDgg0  (3.10) 


by  (3.8). 


-v 
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Along  this  negative  ohareoteri  rtlo ,  therofore 


1 

I 

: 

f 


2  4j“>  ■  2  lo«  J**(*o  io)'  “  "  °'  ■  ^  ^ 8 ♦  higher  order  taros  (3.11) 
'So  *o  2o 

by  differentiation  with  reapeot  to  Rg ,  keeping  oonetant, 

A  firat  approximation  to  the  intermediate  integral  relation  (2.23)  la 
nfi  m  (a  -  fl)6x/r 
or  dJJj)  ■  ed£ 

whence,  by  integration,  and  using  the  boundary  oondltlon  (3.5)  on  fi  at  the 
leading  shook, 


-  -  (*°j>Bq-  (log  g.  -  log  R.) .  (3.12) 

S-A 

Thu a, 

fin-  —  0  3-fiL  a  higher  order  terns  (3.15) 

£' 
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i 


With  thaM  results  the  cycle  may  b«  repeated  to  a  higher  order  of  aoourao]r» 
Thus  a  better  approximation  to  the  relation  (2.21 )  is 


da  +  (a  -  fi  -  (1  -  X*)s*]  S  .  o 

V 


d(“r)  ■  [/»♦  (l  -  X*)a*] 


giving 


n-i'U'-VU- jpftSBpSie. .  x’-yvy 


girlrg 


. .  LSJa .  !SJ5£  h|..M,  *i?WA .  o .  tf&pl 

♦  nigher  order  terms .  (3.1V) 

A  better  apprcxlcetioii  to  the  equation  of  the  posltire  eharaoteristlo  ehioh 
swats  the  leading  shook  at-0  la 


dr(l  -  a  ♦  (l  -  2X*  )fi  *  a*]  -  dt  »  0, 


giving 


*  [, .  isu. .  .  ■?>,..  go)<  h 

.fcfflWl...  0. 

E*  -J 


ehioh  integrates  to  give 


£-t  -  2  •'»„  log£4  2(1  -  X*)(^  gc),[-8  *«  j  108  *  - 

-  2(3  -  2X*)«0',R0  log  £  -  (9  -  VX*)  ^■*flxSq)-- 
■•o-a'iolologSo  -  2(5’-  2’f)*i*S<,  log  So  -  <11  -  «*,)<fe'1S0. 
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Sinoe  fi,  a*  given  by  (3.13)  1»  only  of  the  aame  order  as  the  lower  trder  term# 
of  a  siren  by  (3 *14)  ,  It  Is  not  necessary  to  Improve  the  value  of  fi  In  the  —me 
tray.  A  better  approximation  to  the  equation  (3,9)  for  the  negative 
chnraoteristle  which  lnterseota  the  leading  shock  at  A  la,  from  (2.22), 

dr  [1  ♦  (1  -  2\*)«]  ♦  dt  ■  0. 

Using  the  results  (3,11*.)  and  (3*1 0 1  this  leads  on  Integration  to  the  equation 


£♦1- 


(l  -  2\*)(Cp  Rq)»  (log  Rp  -  log  r) 


2  h  m  V 


The  required  solutions  are,  therefore,  given  by  the  four  relations  (3*13)  to 
(3*16),  but  it  still  remains  to  examine  what  restrictions,  if  any,  are  imposed 
on  the,  as  yet,  unknown  function  *(R) ,  with  the  help  of  the  end -condition, 
«'(R)  •  0  as  £  *  «*,  which  has  not  yet  been  used. 


4*  Bxoess  energy  and 


lse  per  unit  area 


from  the  solution  for  the  deoaying  spherical  shock  given  in  the  last 
seotion,  it  is  possible  to  oaloulate  the  impulse  per  unit  area  of  the  flow 
behind  it,  at  least  as  far  behind  it  as  the  flow  is  free  from  further  shocks. 
Thus  (see  Figure  l),  the  Impulse  per  unit  area  at  a  distance  R^  from  the 
oentre  of  the  explosion,  from  the  time  TA  at  which  the  shook  reaches  this 
distance  up  to  the  time  t D  when  the  positive  oharsoterlstio  through  0  reaches 
this  distanoe,  is  given  by, 


I  ■  /  (p  -  P0)  4t 


or  pJT  *  i  *  /  (fi  *  1)  At  .  (4.1) 

4 

Now  since  D  is  on  the  positive  characteristic  which  meets  the  leading  ehook  at 
0,  by  (3*15)  or  (3.8) 

h.  "  "  2  *o  80  ^8  5»  ■  •»  -  2  *0  So  K>8  5o 

♦  higher  order  tense  (4*2) 

and,  therefore,  differentiating  with  respect  to  R^ ,  keeping  constant, 
dt-  (•„'*L,*  log  L,)' 

~  ■  - rrr - ”  2(  V  Es  ♦  higher  order  terms.  (4*3) 

*0  &  ^ 

The  relation  (4*1)  nay  now  be  integrated  aa  follows! 


--ijdt.ty  (2.8) 
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This  is  a  result  of  great  Importance,  for  It  allows,  at  last,  some 
consideration  to  be  given  to  the  possible  forms  of  the  unknown  function  «(R), 
by  examining  wiiat  happens  to  the  impulse  per  unit  area  I ,  given  ty  (4*4) ,  as 
the  point  C  on  the  leading  shook  tends  to  infinity.  ~ 

First,  in  order  to  proceed  at  all,  it  is  necessary  to  neglect  the 
possibility  of  having  an  infinite  number  of  shocks,  behind  the  leading  shook, 
which  overtake  it.  In  the  absenoe  of  such  an  infinite  array  of  aeoondaxy 
shocks,  it  follows  that,  for  all  values  of  R  greater  than  soma  R^,  there  Is  a 
finite  region  of  shook-free  air-flow  behind  the  loading  shook.  It  may  be 
remarked  here,  in  support  of  this  assumption,  that  observations  do  not  Show  ary 
secondary  shocks  at  ril  which  overtake  the  leading  shook, 

There  are  then  three  possibilities  for  the  form  of  the  unknown  function 

«(R). 

First,  if  «(R)  were  a  function  such  that 

as  R  *  m,  R*  log  R  •  • 

then,  by  (4.4) ,  the  it.  lae  per  unit  area  in  the  air-flow  behind  the  leading 
shock  would  also  ten!  to  infinity  at  some  finite  subsequent  time,  at  all  large 
distances  R^.  For  (see  Figure  l)  ary  further  shook  whioh  Intervened  at  a 
time  tjj,  corresponding  to  a  positive  oharaoteristie  meeting  the  shook  at  a 
finite  point  C,  oust  necessarily  overtake  the  leading  shook  earlier  than  the 
point  0,  and  this  la  a  possibility  whioh  has  already  been  neglected.  The 
possibility  that  »'*  R*  log  R  *  ",  whioh  lead*  to  infinite  positive  legalises 
per  unit  area  at  finite  dlstanoee,  la  also  neglected. 
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Swonl,  if  *(R)  wera  a  function  auoh  that 

M  B  *  “»  a*  log  s  •  0, 

then  the  function 

•i'So’toBSo-loR  *g) 

would  be  eero  for  ^  ^ 

positive  for  So  > 
end  would  •  0  ee  Be  *  "• 

It  would  therefore  have  a  maximum  at  aone  finite  value  g^.  However,  by 
(4.3) ,  at  the  point  Da  corresponding  to  0,, 

tk  m  fa*sfo««Bo  0 

«0  »i  Bo 

to  the  order  of  equation  (4.3) .  ainoe  tj1  g<J(log  B,  -  log  haa  a  marlw 

value  at  D„,  and  beyond  Do,  dtjj/dJ^  would  become  negative.  Thie  possibility 
la  therefore  alao  neglected. 

There  recwina,  then,  only  the  third  possibility  that  t(R)  la  a  function 
auoh  that, 

aa  B  *  «'* B*  1°8  B  ■*  *  finite  non- aero  quantity,  aay  k*/4.  (4.5) 

Then,  aa  B  *  * » 


- * , 

2  E  *>«*  R 

(4.6) 

«  *  k  iog|i  B  ♦  K  • 

(4.7) 

end  it  is  seen  that  this  possibility  ia  consistent  with  the  end-condition  that 
,1  •  0  u  R  *  •.  With  this  possibility,  the  expression  (44,)  for  tha  lapulaa 

per  unit  eras  becomes,  at  R  •  • 

*0 

(l  ♦  X*)k* 

d!A 

(4.8) 

end  thus  as  g^  *  •  , 

I  g*  «  (1  ♦  X*)k*/2  . 

(4.9) 

Turning  now  to  eonalderationa  of  t)w  energy  In  the  air-flow  behind  the 
leading  shook,  the  energy  transported  enroll  the  surface  of  a  given  sphere 
r  ■  from  the  tine  (Figure  l)  at  which  the  shook  reaches  this  sphere  up  to 
the  tine  tp  when  the  positive  eharaoterlstlo  through  0  re  in  he ■  thie  sphere, 
oonslats  of  the  energy  of  the  mast  transported  serosa  this  sphere,  together 
with  the  work  done  st  the  rurface  of  the  sphere  by  fluid  Inside  the  sphere  at 
any  time  on  fluid  out aide  the  sphere  st  any  time.  It  ia  thus  given  by 
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Again,  the  integration  may  be  performed  by  similar  me the -da  to  thoae  used 
previously,  after  deriving  an  expression  for  dj/dJ^  along  A'B  (Figure  1),  and 
the  result  la 


S,  -  4r'^xt  h  *0*  Eo'doa  K0  -  log  g*) 

a  higher  order  terns  (4..  15) 


in  agreement  with  (4. 11). 

The  arguments  applied  to  the  positive  lcpulse  per  unit  area  I,  (4.4.),  in 
order  to  determine  the  moat  probable  fora  of  the  function  «(R),oould,  of  ooursa, 
have  been  applied  equally  well  to  the  energy  B,. 

5.  The  limiting  positive  oharaoterletlo  terminating  the  abode -decay  chase  of 
the  air-flow  behind  a  blast  shock 

In  this  seotion, scene  further  consequences  of  the  reatriotions,  (4.5), 

(4.6),  (4*7) »  plaoad  on  the  function  *(R)  in  the  previous  seotion,  are 
examined. 


Under  these  reatriotions,  the  equation  (5.15)  or  (5.8)  of  the  positive 
characteristic  irliioh  overtakes  the  lending  shook  at  0  reduo e a  to 


r  -  t  -  k  log  j/log^  Rq  ♦  higher  order  terms 


■  K  ♦  0 


(5.1) 


In  the  limit  as  R„  «  **  then,  this  gives  (of.  3>15),  for  the  limiting 
positive  characteristic  which  Just  fails  to  overtake  the  leading  shook  in 
finite  distance,  the  equation 

r  -  t  ♦  k*(l  -  X*yar  -  K.  (5.2) 

This  may  be  eonpared  with  the  result  obtained  by  Whitham  [3]  for  the 
limiting  oharaeterlatla,  namely,  in  tlva  present  rotation, 

r  -  t  ♦  0  (log  3/r)  ■  K .  (5.2a) 

The  method  used  here  appears  to  show  that  the  term  of  order  log  j/r 
envisaged  by  17hitnam  has,  in  fact,  aero  coefficient,  and  that  the  next 
surviving  term  is  only  of  order  yr,  and  is  determinable  by  the  method  used 
here. 


This  limiting  characteristic  (W  in  Figure  1 )  divides  the  whole  fluid 
flow  behind  t>«  leading  shock  into  two  parts,  of  important  significance.  For 
the  flow  behind  the  limiting  characteristic  VW  cannot  exert  ary  influence  at 
all  on  the  ahook  decay  process,  otherwise  than  through  the  mediun  of  a  shook 
overtaking  VST.  But  ary  shock  nhieh  overtakes  VW  must  overtake  the  leading 
ahook  in  finite  distance,  since  Its  slope  rust  be  greater  at  any  point,  than 
the  slope  of  the  positive  oharsoteristio  there.  The  possibility  of  an 
infinite  number  of  such  shocks  has  been  ruled  out,  and  attention  is  being 
confined  to  distances  greater  than  that  at  which  the  last  of  auoh  secondary 
■hooka,  if  any,  overtakes  the  leading  ahook, 

Henoe  no  flow  behind  the  limiting  oharsoteristio  VW  can  influence  the 
shock  decay  process,  and  the  shook  deoay  process  ia  therefore  controlled 
entirely  by  the  flow  between  the  leading  shock  AC  anl  the  limiting  positive 
characteristic  VW.  This  part  of  the  flow  will  be  called  the  shock -deoay 
phase  of  the  blast  wave. 

The  equation  (5.1 )  of  the  leading  ahook  la,  by  (4.7), 

£  ■  T  ♦  k  log^  I  tit  «  higher  order  terma . 


14. 


(5.5) 


The  peek  overpressure  imnediately  behind  the  leading  shook  is,  by  (4.6) 
end  (J.4). 

P  -  1  a  KCI-*. ^  1  +  higher  order  terns.  (5.4) 

£  login 

The  duration  of  the  shock -decay  phase  at  r  ■  RA  (AZ  in  Figure  1 )  is  given 


Ia  "  h  "2d 


-  K  ♦  k*d  -  AV^]  - 


.  '.toe  T,  is  on  the  Uniting  ehexeoteristio  (5,2). 

Thus,  by  (5.3)* 

a  k  log^  ♦  higher  order  terne.  (3.5) 

Behind  the  leading  blast  shook,  the  partiole  valooity  u  ia,  by  (2,13) •  (3.13) » 
(3.14)  end  (4.6). 

tj  a  (l  -  X* ) Vx  logi  Rq  a  higher  order  terms  (3.6) 

and  the  overpressure  (j>  -  1)  is,  by  (2.8),  (2.14) >  (3.13).  (3*14)  and  (4.6), 

£  -  1  a  k(l  ♦  X*)/r  logi^o  ♦  higher  order  terns.  (5.7) 

At  the  limiting  characteristic,  (3.14)  gives,  as  ^ 

a  ■  — -J  ♦  possible  higher  order  tense ,  (5.8) 

and  (3.13)  given,  ae  JL  «  - , 


8m-  ——  +  possible  higher  order  terns.  (5-9) 

4r* 

This,  at  the  limiting  characteristic,  the  overpressure, 

p  -  1  a  (1  ♦  X*)(«  ♦  fi)  ♦  higher  order  terms, 

is  of  lower  order  t.ian  l/r>,  and  the  partiole  velocity  is 

u  a  (l-X*)(«-/»)  a  *  (^  ~  X 1  .  possible  higher  order  terms.  (5»10) 
2r* 

Since,  at  the  limiting  characteristic ,  the  overpressure  ie  not  much 
different  from  aero,  the  si  100k -decay  phase  ee  defined  here  ie  relatively  not 
suoh  different  from  the  pooitive  overpressure  phase  of  the  flow.  However,  if, 
in  feet,  the  overpressure  at  the  rear  of  the  ahock-deoay  phase  is  still 
greater  than  aero ,  it  is  possible  for  e  secondary  shock,  behind  the  limiting 
characteristic ,  anu  never  overtaking  it,  to  enter  the  positive  overpressure 
phase,  increase  the  overpressure  end  thue  prolong  the  positive  overpressure 
phase.  In  explosions  from  conventional  explosives,  such  a  seoondaxy  shook  la 
observed  (5,  6]  and  ita  inception  and  development  have  bean  investigated 
(7,  8). 


t  ,v  % 


To  coi<n>lete  tlie  results  for  the  shock-decay  phase  of  the  motion  It  may  be 
repeated  that  the  Impulse  per  unit  area  In  the  shock -decay  phase  (v.8)  la 


X  ■  . +  higher  order  terms  | 

2R 


(5.11) 


and  the  excess  energy  over  ataospherio  In  the  shock  ml  eo ay  phase  (4.1 6)  Is 

£«  ■  —  *  ■)  ^  «.  higher  order  terms.  (5*12) 

6.  k  solution  to  higher  order  of  aoouracy 

The  considerations  of  section  4  led  to  the  oonelnsion  that  the  only 
restriction  on  the  function  e(R)  for  a  spherically  deoaying  shook  is 


8  *  •> 


«* 


2  a  log*  £ 
k  log^  B  ♦  K. 


There  seems  therefore  no  possibility  of  determining  the  funotion  «  more 
precisely  without  appeal  to  the  initial  conditions  of  the  explosion,  or  to 
boundary  conditions  within  the  leading  spherical  shook,  or  without  making  soma 
hypothesis. 

The  boundary  conditions  at  the  leading  shock  (5.5)  are  axpresnlble  aa 
power  series  in  *'  when  * '  la  mall.  On  Idle  other  hand,  the  funduoental 
length  L  is  so  far  arbitrary,  and  if,  for  instance, a  different  length  L,  were 
used,  and  R ,  ia  written  for  R^L,, 

a.  •  {£) .  [w  i .  *  {£}]* 

.^s[,.^l-i^d . .]  «.,) 

ao  that  an  expansion  in  powers  of  l/l°g  &  would  be  introduced  into  the 
expressions  for  «  and  ,  hy  »  change  in  the  fundamental  length  L. 

On  this  basis,  the  hypothesis  is  made  that,  for  spherics!  explosions,  t 
may  be  represented  by  a  double  series  of  the  form 


e  ■  k  log^  8  “  B  B  “  1  1°*  *  '*  8  ♦ 
♦  i  t .(log*  B) 

♦  *  .(*>«*  8) 


(6.2) 


in  which  the  origin  of  tirar  has  been  oho sen  sc  that  K  n  0. 

It  follows  from  tliia  hypot’iesla  that,  at  large  distanoea  R,  an  attaint  to 
derive  a  higher  order  approximation  than  from  the  first  term  only  of  «,  by 
Including  a  finito  number,  say  two,  of  higher  order  terras,  must  be  directed  not 
towards  including  terms  of  the  order  c'  ani  t'1  times  the  lending  term  but 
rattier  towards  an  approximation  of  the  form 
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i 


t  ■ 


(6.3) 


k  log*  H  -  m  log”*  R  -  1  log”* R  ♦ 


Now  L,  the  fundamental  length,  la  arbitraxy,  and  it  follows  from  (6.1) 
that  L  nay  be  oho  sen  definitely  so  that  the  ooeffieient  m  ■  0.  With  this 
definite  value  of  L  then,  the  approximation  beoomes. 


•  «  k  log*  a 


fi  —VI 

L  k  log*  Rj 


and  all  further  remits  given  are  baaed  on  the  above  hypothesis,  and  with  this 
definite  ohoioe  of  L.  The  expression  (6.4)  oontains,  then,  three  constants, 

L,  1,  and  k  which,  if  the  hypothesis  is  Justified,  should  be  determinable  by 
oocparison  with  experimental  observations. 

With  the  form  (6.4) ,  and  neglecting  terms  of  higher  order  than  l/log*  R 
times  the  leading  term  as  well  as  all  terms  of  order  *'  or  higher,  compared 
with  the  leading  term,  the  results  for  the  shock-decay  phase  of  the  blast  motion 
at  a  distance  R  from  the  centre  of  explosion  aret 


e#  ■ 


r,  *  . 

R  log*  R  L  k  log*  Rj 


Shook  peak  overpressure  (3*4) « 

p-1.  IsLllpc 1+_3L 
“  R  log*  R  l  k  log* 

General*  overpressure  in  the  shock -deoay  phase  (3*7) , 


k  log*  £J  * 


p  1  ■  fi  ♦ _ 21 

£2oS*Sol  kl°6 


8*»j: 


Duration  of  shock-decay  phase  (3*5), 


r  ■  k  -log*  R  fi - i—— ' )  . 

”1  R  log*  Rj  * 


Shook -deday  insulae  per  unit  area  (5.11), 

I  ■  k*(l  +  X*)/®  with  no  term  in  l/R  log*  R  ;  (6.9) 

Equation  of  positive  characteristic  which  overtakes  leading  shock  at  0,  (3.15); 


r  -  t  - 


k  log  r 


log*  So 


['•s*d 


(6.10) 


Relation  between  peak  overpressure,  shock -decay  impulse  per  unit  area,  and 
duration  of  shook^eoay  phase,  (6.6),  (6.8)  and  (6.9), 


(P  -  l)r  a  21  £l  + 


21 
k  log* 


(6.11) 


7.  The  form  of  the 


ee  in  a  spherical  blast  wave 


The  results  givn  in  the  previous  seotion  6  enable  an  expression  to  be 
obtained,  for  the  first  time  so  far  as  the  author  is  aware,  for  the  ultlmte 
pressure -time  relation  behind  a  spherical  blast  wave,  at  sufficiently  large 
distances  from  the  centre  of  the  explosion. 


For,  from  (6.6),  the  peak  overpressure  at  the  shook,  when  it  is  at 
distance  R^ ,  is 


P*1  ■  ft  ♦  — 21 - *1 

k  log*  R»  J 


and  subsequently,  behind  the  shook,  at  the  same  distance  R^  from  the  centre  of 
the  explosion,  the  overpressure  is  generally,  by  (6.7) » 


£  *  1 


.—a 

Sx  1°8*  So  ^  k  log 


id 


at  a  time  tjj  (see  Figure  l)  at  whioh  the  positive  oharaoteristio  overtaking  the 
leading  shook  at  0  reaches  the  distance  R^. 

Since  D  is  on  the  positive  oharaoteristio  whioh  overtakes  the  leading 
shook  at  0,  by  (6.10), 

(7,5> 

whilst  if  r  is  the  shock-decay  duration  AZ  (Figure  l),  by  (6.8), 


r  -  k  log*  R.  fl - 1—)  , 

-*1  klog«Rj 


and  the  relation  between  R^  and  at  the  leading  shock  is,  by  (3»l)  and  (6. if) , 


S*  -  ^  «  k  log* 


From  (7.3)  and  (7*5)*  by  subtraotion, 

^  -  k  log*  R^  1  -  ---•  ^ - 

l  WL  1°8,/,Ra 


- (1  + 


3  log»'»  R.  if  log* 


log"  1 R  log 


*1 '21k, 

S'»R  * 

-a  ./ 


Thus,  from  (7>if), 


(.  -  J  log*  L  k  log* 


Bo  k  log  R^  log  Rq  k  log1 


’8*  Ra]’ 


whereas,  from  (7.2)  and  (7.1 ), 

p  -  1  log*  ga 

p^t  *  — nr 


.'•og a  s. 


.i+—4--_ i i-] 

k  log*  ^  klog»  R^j 
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Thu*,  from  (7.7)  aid  (7.6), 


BH' 

to  th«  present  order  of  approximation. 

If,  In  tho  present  oontext,  time  la  measured  from  the  instant  the  shook  reaches 
the  radi.ua  R^,  (7*9)  nay  he  written  siqply, 

‘  (wo) 


This  may  he  oompared  with  the  empirical  form  first  suggested  by  Friedlander 
[9],  namely. 


H*  ■  MK 


or  with  Tariationa  on  the  Friedlander  form  [of.  5,  6]  such  as 


S  ■  Nv,v- 


p  -  p, 


*  -p, 


(7.11) 


(7.12) 


(7.13) 


In  the  relations  (7.1l)»  (7*12)  and  (7.13),  o,  and  ea  are  oonstants,  and  T  is 
clie  duration  of  the  positive  overpressure  phase. 

To  obviate  the  experimental  difficulties  associated  with  the  measurement 
of  quantities  connected  with  the  entire  positive  phase  or  the  entire  chock- 
decay  phase  as  introduced  in  section  5#  it  is  oonveni  to  define  a  general 
time-interval  Tf,  after  the  arrival  of  the  leading  blast  shock  at  a  given 
radius  ,  as  tho  elapsed  time  at  which  the  overpressure  has  fallen  to  a 
fraction  f  of  its  peak  value.  The  positive  characteristic  through  the 
particular  point  D<  (Figure  l)  corresponding  to  (R^,  ♦  rf)  overtakes  the 

leading  shock  at  the  particular  point  0,  satisfying,  by  (7.8), 


log^R. 

f  ■  -.rn. 


whilst,  from  (7*6), 


k  log* 


^  r, .  *4^-]  -  — 2 _ r 

l  ^SoJ  Wfc|;  iog,/i  Hg(  W'* &>,  y 


i  Iclog  R^  klog'R^j’ 


(7.14) 


From  (7.14)  and  (7.15),  to  the  present  order  of  approximation. 


(7.15) 


(7.16) 
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Of  particular  interest  is  the  epeoial  case  whan  the .  fraction  f  has  the 
value  fx  ■  (Vj  -  i)/2.  In  this  particular  oaae 


(1  ♦  2fa)  *  .  3 


(7.17) 


k(3  -  ^3) 


log*  &»  <  1  -  ■■  ■ 

k  log  H.  j 


(7.18) 


so  that,  from  (7*4) » 


■f! .  main _ 2lji 

I  2  l1  klog»*J 


(7.19) 


and  (6.1 1)  may  he  written 


(P  -  1  )lx  -  (3  -  *>)  I. 


(7.») 


The  relation  (7.20 )  appears  to  be  the  moat  useful  relation  of  its  kind,  for 
further  analysis  gives,  in  general, 


.  kin ±j£HLi£!iri JS£L-1 , 

L  kkg^J' 


(7.21) 


as  a  generalisation  of  (5.11 ),  and  thus 


(£-0  lt  -  7-— -r  fi  ♦.2-10-Z-f).*.1l 
(i  ♦  t)  l  k  iog“  ^  J 


(7.22) 


as  a  generalisation  of  (6.11).  It  is  olear  from  (7.22)  that  the  speolal  oase 
f  ■  f jj  has  no  particular  merit  so  far  as  this  relation  is  oonoemed,  and  that 

there  is  no  other  speoial  value  of  f  (apart  from  the  trivial  exception  f  ■  l) 
which  will  reduce  (7.22)  to  the  form 


(S-^Jf 


constant. 


Finally,  the  generalisation  of  (5*12)  is  obtained  as, 

$,f  -  Ai  -  f)  B,  [i  -  •  (7.23) 

son  with  experimental  observations 

In  order  to  facilitate  oorparison  with  experimental  observations,  it  is 
desirable  to  derive  some  further  results  of  the  theory  in  terms  of 
physical  quantities  rather  than  the  reduced  non-dimensional  quantities  which 
have  been  uaed  in  the  development  of  the  theory. 

Thus,  equation  (6.9)  may  be  written 

kV  .  .  (a  4 


*  X«)  * 


end  equation  (7.20)  may  ‘>e  written 


k*L*  .  -.0  V*»o  Ofr  ♦  1) 

(1  ♦  X*)  /3 
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in  which  (P  -  l)  ■  (P  -  Pq)/p0  ia  the  peak  overpressure  measured  in 
atmospheres  and  is  left  in  the  reduced  notation  for  sinplioity. 

From  (6*6)  and  (6.8) 

(P  -  1)  r*  R  ■  k4(l  ♦  X*)(logR), 
and  this  may  now  be  written 


kV  (log  R  -  log  L)  ■ 


(P  -  1)  rs  R  a* 


(1  ♦X1) 


a-. 


Similarly,  from  (6. 8)  and  (V.18) 


r> 

T 


(8.3) 


and  this  may  now  be  written 


kV  (log  R  -  log  L) 


**  *g(3  -  V$) 


2r„ 


Eliminating  kL  between  (8.3)  and  (8.4)  gives  the  result 

log  L  ■  log  R  -  3„(2--  ll  . 

8  8  2  8^  -1)  rj 


(B.4) 


(8.5) 


Finally,  when  these  five  relations  have  been  found  to  yield  consistent  values 
of  k  and  L,  by  comparison  with  experimental  observations ,  the  following  three 
relations  may  be  used  to  determine  1,  and  thus  further  assess  the  consistency 
of  the  agreement  between  theory  and  observations. 

From  (6.6) 


From  (6.8) 

1  «  k  log*  R  -  r  log*'*  R,  (8.7) 

and  from  (7.18) 


(P  -  l)  R  log*'*  R  k  log*  R 
1  "  3(1  ♦  X*)  ‘  3  “* 


(8.6) 


1 


3  3V3 


tm  log"*  R. 


(8.8) 


Potter  and  Jarvis  (5,  6]  have  made  observations  of  explosions  from  bare 
spherical  charges  of  TUT  and  RDVTOT  6O/4D  which  are  sufficiently  comprehen¬ 
sive  for  a  comparison,  on  the  above  lines,  with  the  theory  developed  here. 
Their  results  are  given  in  terms  of  the  oonxoon  method  of  scaling,  using  a 
factor  TT3  (W  ■  charge  mass  in  pounds).  Here,  W  is  put  equal  to  one  pound  in 
their  jresults,  so  that,  in  faot,  their  observations  are  aoaled  to  one  pound  of 
explosive  at  sea-level,  and  the  comparison  is  made  on  this  basis.  On  this 
scale,  their  observations  extend  from  about  4  ft.  to  50  ft.  from  the  centre  of 
the  explosion,  and  their  peak  overpressure  measurements  were  fitted,  over  the 
oooplete  range,  by  the  method  of  least  squares,  to  the  form 


P  -  1  >  v*R  ♦  Vk*  ♦  <^R». 

Potter  and  Jarvis  have,  of  oourae,  analysed  their  records  in  terns  of  the 
usual  positive  overpressure  phase.  In  the  oaae  of  RDVtNT  they  observed  a 
secondary  shook  which  entered  the  positive  phase  at  about  15  ft.  from  the 
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centre  of  the  explosion,  and  at  greater  distances  their  quoted  values  of 
positive  duration  inolude  a  email  extrapolation  to  give  an  estimate  of  what  the 
positive  duration  woull  have  been  in  the  absence  of  a  second  shook.  Their 
values  of  positive  impulse  per  unit  area,  however,  inolude  only  the  inpulse 
in  the  primary  wave.  The  comparison  made  here  rests  entirely  on  the 
assumption  that  these  results  of  Potter  and  Jarvis  may  be  taken  as  a  good 
approximation  to  values  for  tlie  shock -decay  phase  introduced  in  the  present 
theoiy. 

It  may  seem  surprising  that  the  results  of  Potter  and  Jarvis  are 
sufficiently  comprehensive  to  include  measurements  of  (7«l8).  In  fact. 

Potter  anl  Jarvis  have  measured,  on  their  records,  the  time  ,  (7*1 6) ,  for 
the  fraction  f  ■  l/e.  By  a  remarkable  numerical  coincidence, 

J_  .  _L_  .  /j  +  1  «  2.73205, 
f x  *3-1 

as  compared  with  e  ■  2.71828,  or 

(1  +  §)*  a  3.0129, 
as  oonpared  with  (1  +  2fx)*  ■  3» 

The  values  of  r t  /e  ,(whioh  they  name  the  'decay-constant' ) ,  given  by  Potter  and 
Jarvis,  have  therefore  been  used  as  values  of  r„.  The  author,  of  course, 
disagrees  violently  with  their  assertion  that  'the  ohoioe  of  the  factor  l/®  !■ 
quite  arbitrary,  and  0.5  or  0.25  would  have  been  almost  aa  good'! 

Table  (8.1 )  shows  the  results  for  one  pound  of  RD^TOT  60/40  at  so  a -level. 
Columns  (6)  anl  (7)  show  the  two  values  for  k*L*  given  by  the  relations  (8.1) 
aid  (8.2).  Column  (8)  ahows  the  value  of  L  given  by  the  relation  (8.5). 

Figure  2  shows  the  plots  of  the  relations  (8.3)  and  (8.4),  on  semilcgarlthmio 
graph  paper,  anl  the  straight  lines  corresponding  to  the  values  k  ■  0.718, 

L  ■  3.8  ft.  are  also  drawn.  Finally,  using  these  values  of  k  and  L,  the  three 
values  of  1  given  by  the  relations  (8.8),  (8.7)  and  (8.8)  are  given  in 
oolumna  (9),  (10)  and  (ll)  of  Table  8.1,  being  denoted  respectively  by  l(P), 
l(r)  and  l(rx). 

Table  (8.2)  anl  Figure  3  give  the  corresponding  comparison  with  the 
results  for  one  pound  of  TOT  at  sea-level,  and,  for  TOT,  the  values  k  ■  0.8025, 
L  ■  4*5  feet  have  been  chosen.  At  the  greatest  distances  from  the  oentre  of 
the  explosion,  the  observations  naturally  become  less  reliable,  and  the  least 
squares  fit  of  the  mean  curve  by  Potter  and  Jarvis,  being  chosen  as  an 
optimum  over  the  whole  range  of  observations,  is  probably  not  so  good  at  the 
greatest  distances,  as  a  best  fit  over  a  shorter  range  of  larger  distances. 

It  is  not  surprising,  therefore,  that  the  good  agreement  shown  in  Tables  8.1 
anl  8.2  aid  Figures  2  anl  3  up  to  distanoes  of  about  35  feet,  ia  not  so  well 
maintained  at  4D  and  50  feet.  It  appears  that  the  fitted  moan  values  of 
(P  -  l)  and  I  are  the  least  reliable,  and  amended  values  of  these  quantities 
are  shown  in  brackets  in  Tables  8.1  and  8.2  and  Figures  2  and  3,  to  give  sane 
idea  of  the  changes  in  these  quantities  which  would  give  agreement  with  the 
theory  given  here.  Those  changes  are  quite  small. 

Finally,  as  an  alternative  to  (8.2),  the  relation  (7.20)  may  be  written 
in  the  form, 

(P  -  P0>  r*  -  (3  -  ^)  I.  (8*9) 

anl  thus  constitutes  an  overall  check  on  the  hypothesis  of  seotion  7»  when 
tested  against  experimental  observations. 

The  values  of  (P  -  P0)rj/(3  “  ^3)  I  are  given  in  column  12  of  Tables  8.1 
end  8.2.  It  is  seen  that,  in  both  oases,  as  R  increases,  this  quantity  tends 
rapidly  and  monotonicaUy  to  unity. 
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9.  Remarks 

(a)  It  needs  to  be  streased  that  the  hypo the sis  made  in  section  6  and  the 
consequent  results  given  in  seotions  7  and  8  are  not  necessarily  sound,  and 
other  similar  hypotheses  need  to  be  investigated  before  any  firm  conclusions 
can  be  drawn.  From  the  results  given  in  section  8,  however,  it  is  seen  that 
this  hypothesis  appears  to  be  in  reasonable  acoord  with  experimental 
observations,  particularly  in  the  sense  that  the  observations  used  are  far 
more  extensive  than  the  minimum  requirements  for  determining  only  three 
constants  L,  k,  1,  and  yield  consistent  values  whenever  the  constants  are 
determined  by  alternative  methods  from  different  types  of  observation.  On 
the  other  hard.,  the  solution  which  is  being  compared  with  experiment  is  an 
asymptotic  one,  and  the  vital  determination  of  the  oonstanta  depends 
therefore  on  the  measurements  at  large  distances,  where  the  observations 

are  naturally  at  their  weakest.  For  this  reason  the  numerical  oonstanta 
chosen  to  fit  the  observations,  in  section  8,  should  not  be  regarded  as 
determined  to  any  high  order  of  accuracy. 

(b)  It  is  of  interest  to  compare  the  fit  with  observations  of  the  higher 
order  asymptotic  solution  proposed  here  with  the  corresponding  fit  of  the 
•aero  order'  solution  involving  only  the  leading  terms.  Kirkwood  and 
Brinkley  [2  et  al  ] ,  who  gave  a  solution  to  the  whole  development  of  a 
spherical  explosion  by  making  a  physical  assumption,  obtained  the  aaymptotio 
solution  for  the  peak  overpressure  from  the  limiting  form  of  their 
differential  equations,  as  the  leading  term  only  of  the  solution  given  here, 
i.e.  in  the  present  notation 


P.i  .  Ml.-JiV. 

B,lag*R, 

in  terns  of  a  fundamental  length  Li ,  where  R,  ■  l/L, . 


(9.1) 


The  corresponding  solution  for  other  quantities  to  the  same  order  of  accuracy, 
by  the  methods  given  here,  would  thus  be 


k?(l  ♦  X«) 


®i 


r  ■  ki  log?  R, 


(9.2) 

(9.3) 


and,  in  fitting  this  solution  to  experimental  observations,  only  two  oonstanta, 
k,  an’  Li ,  have  to  be  determined.  Kirkwood  and  Brinkley  do  not  discuss  the 
physi*  .1  significance  of  their  length  L,,  but  fit  the  solution  (9.l),  (9*2), 
(9«3)  asymptotically  to  the  numerical  solution  of  their  equations  closer  in, 
and  thus  derive  definite  values  for  k,  and  Li*  In  faot,  when  written  in  the 
form 


p-p0  (k*L,)(1  ♦  \«) 

po  R  (log  R  -  log  L,)*  * 

(9.4) 

I  «o  (kiL,)«(l  +  X*) 

(9.5) 

P0  2R 

ra0  .  (k,Li)(logR  -  log  L,)*  , 

(9.6) 

it  is  olear  that  this  solution  oan  be  used  to  determine  definito  values  of  k, 
and  L,  by  fitting  to  experimental  results  or  calculations,  and  thus  the  form 
of  the  solution  (9.4)  to  (9.8)  involves  inherently  some  definition  of  a 
particular  Li,  different  from  the  preoise  definition  of  L  used  in  the 
hypothesis  of  4*«otion  6. 
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Since  equation  (9.5)  is  identical  with  (8.1 ),  the  fitting  of  thia  solution  to 
the  same  set  of  observations  should  give  the  same  value  for  to  Li  as  obtained  in 
soot  ion  8  for  kL. 

For  oonparison,  some  values  of  k,L,  obtained  by  Kirkwood  and  Brinkley, 
both  by  choosing  the  initial  data  of  their  solution  to  fit  experimental  data, 
and  by  choosing  initial  data  from  the  thermodynamio  properties  of  the 
explosive,  are  quoted  here. 

Referenoe  2.  1  lb.  Cast  Pentolite 

From  peak  overpressure  solution  k?L*  a  8.73  ft? 

(L,  a  8.09  ft.,  R,«  0.365) 
From  positive  insulae  solution  k*L*  a  8.53 

Reference  2.  1  lb.  Torpex  II 

From  peak  overpressure  solution  k?L*  a  9.21  ftf 

(L,  a  8.32  ft.,  k, . 0.365) 

From  positive  iopulae  solution  k?L*  a  8.65  ft? 


1  lb.  TNT 

From  positive  impulse  solution 
From  positive  overpressure  solution 


kfrf  .  7.97 
k*L*  a  7.96 
(Li  a  8.41  ft., 


k, -0.335) 


Reference  18.  1  lb.  TNT 

From  positive  impulse  solution  a  7.  JO 

From  positive  overpressure  solution  k^L*  a  7.49 

(L,  -8.14  ft.,  k,- 0.336) 

(o)  It  has  frequently  been  suggested  that  a  blast  wave  may  ultimately  tend  to 
the  form  of  an  N-wave  with  a  second  shook  about  equal  in  strength  to  the 
leading  shock.  The  development  given  here  suggests,  on  the  other  hand,  that 
no  conclusions  oan  be  drawn  about  the  air  motion  behind  the  limiting 
oharaoteriatio,  on  the  strength  of  the  governing  equations  and  leading 
boundary  conditions  only.  Experiment  [5*6  ]  shows  that  the  shook  corresponding 
to  the  second  leg  of  the  N  is  consistently  very  much  weaker  than  the  leading 
shock.  In  the  case  of  a  conventional  explosion,  moreover,  the  second  leg  of 
the  N  ia  the  third  abode,  and  the  shock  next  behind  the  leading  shook  ia 
roughly  in  the  middle  of  the  sloping  part  of  the  N.  It  ia  thia  aeooni  shook 
which  [of.  7,8]  originates  at  the  aurfaoe  of  the  explosive  and  ia  reflected 
outwards  again  after  implosion  at  the  oentre.  This  type  of  shock  is  peculiar 
to  conventional  explosives,  and  consequent  upon  the  infinite  pressure  gradient 
immediately  behind  a  spherical  detonation  wave. 


It  is  the  author's  opinion  that  secondary  shocks  are  not  essential  to  the 
ultimate  forms  of  blast  waves,  and  that  their  presence  or  absence  is  dependent 
on  the  nature  of  the  explosion,  or,  more  precisely,  on  the  motion  of  the 
'explosive  products'  which  sot  like  a  piston  in  causing  the  explosion.  As 
remarked  above,  the  second  shock  from  a  conventional  explosive  is  peculiar  to 
thia  type  of  explosive*  In  atomic  explosions  no  secondary  shook  is  observed 
at  all  Li 5,  P&ge  50,  footnote].  The  author  believes  that,  if  an  explosion 
were  oaused  by  a  spherical  piston  which  returned  to  its  initial  position,  a 
true  N-wave  would  ultimately  form,  Just  as,  in  the  corresponding  case  of  a 
finite  body  in  steady  superBonir  flow,  the  bow  and  tail  shocks  ultimately  form 
an  N-wave  far  out  from  the  body.  In  this  latter  oase,  the  shape  of  the  body 
determines  the  'piston  curve',  whioh,  when  the  body  is  finite,  returns  to  its 
initial  position. 


The  author  is  greatly  indebted  to  Dr.  A.  Coonibs  for  extensive  checking 
and  emendation  of  the  menu script ,  and  discussion  of  many  details  in  the 
presentation. 
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